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Abstract 

Evolution of the spatial arrangement of cells in a primary culture of cardiac tissue derived from newborn rats 
was studied experimentally over extended period. It was found that cells attract each other spontaneously to form 
a clustered structure over the timescale of several days. These clusters exhibit spontaneous rhythmic contraction 
and have been confirmed to consist of cardiac muscle cells. Addition of a contraction inhibitor (2,3-butanedione-2- 
monoxime) to the culture medium resulted in the inhibition of both the spontaneous contractions exhibited by the 
cells as well as the formation of clusters. Furthermore, the formation of clusters is suppressed when high concentrations 
of collagen are used for coating the substratum to which the cells adhere. From these experimental observations, it 
was deduced that the cells are mechanically stressed by the tension associated with repeated contractions and that 
this results in the cells becoming compact and attracting each other, finally resulting in the formation of clusters. 
This process can be interpreted as modulation of a cellular network by the activity associated with contraction, which 
could be employed to control cellular networks by modifying the dynamics associated with the contractions in cardiac 
tissue culture. 
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1. Introduction 

Primary cultures of cardiac muscle cells have been 
successfully utilized as experimental models to study 
the complex electrical and mechanical activity of 
the heart. This is because, compared to experiments 
using an intact heart, a variety of methods can be 
employed with relative ease to visualize and mod- 
ify the dynamical activity of cells. Taken together, 
these advantages enabled us to explore the nature 
of the complex dynamics that arise in cell cultures 
[1,2,3,4,5,6]. 
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Previous studies revealed that the spatio- 
temporal patterns in the cellular activity of rhyth- 
mic contraction in cardiac muscle cell cultures un- 
dergo qualitative changes depending on the age of 
the culture. Soen et al. investigated the spontaneous 
contraction of a monolayer of cardiac cell culture 
for over one week, and found that various inter- 
esting temporal patterns, including subharmonic 
and alternant rhythms, arise spontaneously [7]. 
Spontaneous emergence and annihilation of spiral 
waves were also observed in several studies [2,3,8]. 
Kawahara and coworkers reported that formation 
of cellular network decreases the fluctuations in the 
timings of contraction [9,10]. One possible reason 
for this phenomenon is an alternation in the degree 
of gap junction formation which brings a passive 
conductance between cells [11,12]. Interestingly, de- 
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spite having been implicated in affecting the charac- 
teristics of intercellular communication, relatively 
little attention has been directed at the evolution 
of the spatial arrangement between cells. This may 
be due to the current belief that cardiac muscle 
cells, because they differ from self-propelling cells, 
including fibroblasts, do not possess the molecular 
apparatus necessary for facilitating active migra- 
tion. However, even if a single cardiac muscle cell is 
unable to migrate by itself, it may be possible for 
the geometrical arrangement of the cells in a popu- 
lation to change via interactions among the cells. 

To more accurately assess such a possibility, we 
investigated the evolution of the spatial arrange- 
ment of cells in a cardiac cell culture over an ex- 
tended timescale in the present study. We found that 
cardiac muscle cells spontaneously aggregated and 
formed clusters in which cells constituting the same 
cluster periodically contracted in synchrony. The ap- 
plication of an agent that inhibits the contraction 
of cells results in the suppression of cluster forma- 
tion, suggesting that the periodic contraction of in- 
dividual cardiac muscle cells affects cluster forma- 
tion. Furthermore, the degree of cluster formation 
depends on the concentration of extra cellular ma- 
trix to which the cells adhere. These experimental 
findings indicate that the tension generated by peri- 
odic contraction of cardiac muscle cells contributes 
markedly to the driving force underlying cluster for- 
mation. 

The results of the present study suggest that the 
rhythmic contraction of cardiac muscle cells affects 
the spatial arrangement and structure of the cellu- 
lar network. This phenomenon is likely to be con- 
cerned with age-dependent changes in the spatio- 
temporal dynamics of the electric activity of cells 
observed in earlier studies. This new hypothesis, 
contraction-induced modulation of an intercellular 
network, might provide a new class of dynamical 
phenomena in a cardiac tissue culture. 

2. Materials and Methods 

2.1. Preparation of cell culture 

Primary cultures of ventricular cells were pre- 
pared from one-day-old newborn rats according 
to methods described elsewhere [13,14]. Briefly, 
ventricles isolated from the hearts of one-day-old 
newborn rats were minced and were treated with 
collagenase. Isolated cells were collected by cen- 



trifugation and were supplied with plating medium 

(Dulbecco-modified Eagle medium containing 10% 
fetal bovine serum, 1% penicillin/streptomycin and 
0.1 mM 5-bromo-2'-deoxyuridine (BrdU)). In order 
to reduce the population of non-muscle cells, the 
method of selective plating was employed [13,14]. 
Then, the cells were counted and diluted to desired 
densities (from 2 x 10'"' cells/ml to 6 x 10^ cells/ml). 
Because accurate estimation of the density of cells 
at this stage is technically difficult, not the absolute 
vahi(^s but th(^ relative values of initial cell density 
given ill the following sections have a meaning. The 
cells were plated on the collagen-coated petri dishes 
(see below), which was followed by iiicubaticm for 
24 h. Subsequently, the medium was exchanged 
with contraction medium (modified Eagle medium 
with 10% calf serum, 1% penicillin/streptomycin, 
and 0.1 mM BrdU), and the cells were incubated 
under the same conditions as described above. 

2.2. Preparation of collagen- coated culture dish 

Polystyrene petri dishes (BD Biosciences, CA) 
coated with collagen were used for experiments. In 
experiments described in Sections 3.1 to 3.4, dishes 
were coated with collagen by incubating them with 
collagen (Type I, derived from rat tail, WAKO, Os- 
aka) solubilized in acetic acid (0.1 %, pH 3) for sev- 
eral minutes at room temperature. In experiments 
described in Section 3.5, the amount of collagen used 
to coat dishes was controlled as follows: First, vary- 
ing concentrations of collagen solutions (1.2 mg/i 
to 120 mg/£) diluted in phosphate-buffered saline 
(PBS) were applied on petri dishes, which were then 
incubated at 4°C for 7 h. Dishes were washed with 
PBS several times before use. 

2.3. Time-lapse observations 

In order to observe the process of cluster forma- 
tion, time-lapse observation method was adopted. 
After exchanging medium on the first day, a cul- 
ture dish was positioned on an on-stage incubation 
chamber (IN-0NI-F2, Tokai Hit, Shizuoka) and 
examined under an inverted microscope (IX-70, 
Olympus, Tokyo). The culture was maintained un- 
der humidified conditions at 37°C and 5% CO2. 
The cell culture was observed using a lOx objective 
lens under bright-field illumination and images were 
captured every 3 min for up to a week, using a dig- 
ital charge- coupled-device (CCD) camera (C4742- 
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95-12SC, Hamamatsu Photonics, Shizuoka). The 

culture medium was not exchanged during the 
time-lapse observations. 

2.4. Determination of beat-rate distributions 

The spatial distribution of the beat rate of car- 
diac muscle cells was studied as follows. Culture 
dishes were incubated at the same condition as de- 
scribed above. Spontaneous contraction of cardiac 
muscle cells were observed every 24 h using inverted 
phase-contrast microscope (4x magnification) and 
recorded at 30 frames per seconds for 3 minutes 
under the room temperature (25 °C). The sequen- 
tial images (960 /zmx 960 /zm) were subdivided into 
small segments of 40/xmx 40 /J.m (the size roughly 
corresponding one or two cells), for each of which 
average gray value was calculated to yield a set of 
time series. In these time series, clear peaks associ- 
ated with spontaneous contractions of colls in each 
segment occur [7] . Thus, the number of peaks in one 
minute was counted in order to calculate the beat 
rate of cells in each segment. 

2.5. Determination of cellular orientation 

The orientational order of cells was investigated, 
following the procedure of Ref. [15], using bright- 
field time-lapse images. First, an image was subdi- 
vided into segments of 15 /imx 15 /im, and a Fourier 
power spectrum was calculated for each segment. 
The cross-correlation matrix of each Fourier spec- 
trum was calculated, and the principal axes were 
determined. It has been known that the principal 
axis with the smaller eigenvalue represents the mean 
direction of cellular arrangement [15]. The differ- 
ence between the two eigenvalues can be interpreted 
as the degree of orientational order, because it in- 
creases as the anisotropy of the original image in- 
creases. 

2.6. Histochemical analysis 

Cells were stained with a monoclonal antibody, 
5C5, directed against a-sarcomeric actin (A2172, 
SIGMA, Tokyo). Furthermore, Vectashield (Vector 
Laboratories, U. K.) was used to overlay the cul- 
tures to prevent the sample from drying out. Vec- 
tashield also contains 4',6-diamino-2-phenylindole 
(DAPI) which stains the nuclei of the cells. The 
stained cells were observed with an epifluorescent 



inverted microscope (IX- 70, Olympus, Tokyo), and 
images were captured rising a digital CCD cam- 
era (DV887ECS-UVB, Andor Technology, Northern 
Ireland). 

2.7. Analysis of cell distribution 

The spatial distribution of cell nuclei was analyzed 
to assess the spatial distribution of cells. 

First, the distribution of the cellular clusters were 
determined as follows. At the seventh day from the 
onset of incubation, the cells in culture were fixed 
using cold methanol (-20 °C). The phase-contrast 
images of the culture were obtained for 16 x 16 fields 
of view, each of which has the size of 880 /xmx 704 
fim and is 1 mm apart from its neighbors, using 
a motorized mechanical stage (BIOS-105T, Sigma 
Koki, Tokyo). The number of clusters was counted 
manually from these images. 

Second, the distribution of the cells was deter- 
mined at a higher resolution as follows. The fluores- 
cent images of cell nuclei stained with DAPI were 
obtained for at least 79 fields of view (820 /im x 820 
/xm each) for every culture dish. At least 9700 nuclei 
are included in the analysis for each sample. No sig- 
nificant difference among the samples in the same 
condition was observed in the essential trends in the 
results of following analysis. The images were pro- 
cessed and analyzed with image processing software 
ImageJ (a public domain software originally devel- 
oped at the US National Institute of Health). The 
images were converted to binary images by adopting 
a manually determined threshold of gray value. Be- 
cause there is a tendency that the images of neigh- 
boring cell nuclei become indistinguishable in the 
region with high cellular density, i.e., in clusters, the 
watershed algorithm was applied in order to sepa- 
rate them [16] , although the separation was not com- 
plete. This leads to underestimation of the cellular 
density in the vicinity of clusters. 

On the basis of the coordinates of the identified 
nuclei, a radial distribution function for cell nuclei 
was calculated as 



where n{r)dr denotes the average number of nuclei 
falling within a distance of (r, r+dr) from a reference 
nucleus. In order to avoid the boundary effect, in 
each field of view, the reference cells were selected 
from a region distant from a certain distance R from 
the boundary of the image {R = 200 /xm). g{r) has 
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the unit of number density, and is normalized so that 
its value at a sufficiently large r gives the average 
number density of nuclei. 

Third, the distribution of the size of clusters was 
determined as follows. For each field of view, cellular 
nuclei whose center is within 50 /im from the cen- 
ter of another nuclei were considered to be included 
in the same cluster. Then, the number of cells con- 
stituting the same cluster was calculated for all the 
fields of view. 

2.8. Analysis of contractile motion of cells 

In order to characterize cellular contraction, mo- 
tion analysis of fluorescent microbeads attached to 
the surface of the cells was conducted. The culture 
medium was exchanged on the sixth day of cul- 
ture with contraction medium containing 2 x 10~^ 
% (v/v) of fluorescent polystyrene beads (Molecu- 
lar Probes, 1 /xm diameter). After one hour of incu- 
bation, these beads were sedimented in the culture 
dish and became uniformly attached to the cell sur- 
face. The density of the beads was calculated to be 
2.4 X 10^/mm^ (approximately six beads per cell). 
The beads were found to move in phase with the con- 
tractions exhibited by each cardiac muscle cell. The 
motion of the fluorescent beads was captured using 
an intensifled CCD camera (C2400-27, Hamamatsu 
Photonics, Shizuoka). Using Image J software, the 
proportion of beads with centers of mass that had 
trajectories exceeding 3.0 fim were calculated. 

3. Results 

3.1. Formation of clusters 

First, the cells were plated at an initial density of 
4.4 X 10^ cells/mm^. Changes in the arrangement of 
cells over time was recorded using time-lapse obser- 
vation methods; as shown in Fig. 1, the cells exhib- 
ited a variety of dynamic organization structures. 

Initially, (24 h to 44 h from the onset of incuba- 
tion), the distribution of cells was approximately 
uniform, but the number of the cells that adhered 
to the siibstrate increased with time. At this stage, 
a part of the cells started exhibiting spontaneous 
contractions. However, the frequencies and phases 
of periodic contraction were observed to vary among 
cells, and synchronization of periodic contraction is 
restricted to small aggregations of cells (see Fig. 2a) . 
After approximately 54 h of incubation, small clus- 



Fig. 1. Changes in the arrangement of cardiac cells over time. 
The image sequence shows the changes in the cell culture 
observed from the onset of incubation to 134 h. The scale 
bar represents 500 /xm. 

ters of cells became distinct. As aggregation of the 
cells proceeded, they formed a network of intercon- 
nected small clusters (74 h to 94 h). At this stage, 
the cells that formed aggregations were found to 
contract periodically, and the cells belonging to the 
same cluster tended to exhibit a synchronized con- 
tractile motion; however, global synchronization of 
the timings of periodic contraction was not observed 
(Fig. 2b). Aggregation of cells continued, with the 
interconnected network of small clusters forming 
several distinct clusters of cells (after 124 h), and 
acquiring a round shape (134 h). At this stage, 
spontaneous contraction was only observed in the 
large, developed clusters, whereas the cells outside 
clusters did not exhibit any rhythmic contractions. 
The periodic contractions of the large clusters were 
globally synchronized throughout the culture dish 
(Fig. 2c). 

As illustrated in Fig. 3, the distribution of clusters 
in the culture dish was almost uniform, although 
there exist an inhomogeneity caused by the initial 
inhomogeneity of cells. At this initial density of cells 
(4.4 X 10^ cells/mm^) , the number density of clusters 
was 22 ± 5 /mm^ (Mean ± S. D.). 
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Fig. 2. Changes in the beat rate of spontaneous contraction 
through development of clusters. The upper images display 
the phase-contrast microscopic images at a) 24 h, b) 72 h, 
c) 120 h from the onset of incubation. These images were 
captured for the different samples. The scale bar represents 
300 ^m. The lower images represent the beat rates of spon- 
taneous contraction for the same fields of view as the upper 
images. Beat rates (beats per miimtcs) were displayed using 
the gray level indicated on the right of each figure. 
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Fig. 3. Number density of clusters on the seventh day from 
the onset of incubation. The scale bar represents 2 mm. The 
number density (number of clusters per mm^) is represented 
as a gray value indicated at the right of the image. 

Fig. 4a displays the process of cluster evolution 
in the early stage at a higher spatio-temporal res- 
olution than that shown in Fig. 1. As seen in this 
figure, the parts with higher initial density serve as 
nuclei of cluster formation. Fig. 4b also displays the 
images of higher spatio-temporal resolution for the 
later stage of cluster formation. For the same field of 
view, the field of orientational order of cells, which is 
obtained from the local Fourier transform of the im- 
ages in Fig. 4b, is represented in Fig. 4c. As observed 
in this figure, the clusters shrink into the direction 
in which they are oric'nt(;d. It has been known that 
the muscle fibrils, which are the cellular apparatus 
responsible for inducing contraction, tend to align 
along the long axis of the cell and contraction occurs 
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Fig. 4. Process of cluster formation at a higher resolution 
both on space and time for the same experiment as in Fig. 1. 
The scale bar represents 100 /im. a) The early stage of cluster 
formation. The time from the onset of incubation is displayed 
on the each image (see also the supplementary movie), b) The 
later stage of cluster formation (see also the supplementary 
movie), c) Orientation of the cells is represented by solid 
lines for the same field of view as in b). The angle of ea<;h 
line represents the direction of the principal axis with the 
smaller eigenvalue of the cross-correlation matrix of the local 
Fourier transform. The length of each line is proportional 
to the difl'erence between two eigenvalues, which therefore 
represents the degree of orientational order. 



5 



a) 





b) 




























































I 


































- 

















































































































































Fig. 5. a) Fluorescent image of cardiac cell culture on the 
sixth day after the onset of incubation represented in pseu- 
docolor, a-sarcomeric actin is colored red, while the cell nu- 
clei are shown in blue. The scale bar represents 100 fim. b) 
Local density of nuclei calculated for the same field of view 
as shown in a). Density, displayed as a color map, is repre- 
sented as the number of nuclei in each small block (6.6 X 10^ 

in this direction. The observation given in Fig. 4b 
and c thus imply that the contraction of cells is re- 
lated to the process of cluster formation and this 
possibility is discussed below. 

3.2. Identification of cell type 

Given that a primary culture was employed in the 
present study, the culture was contaminated by non- 
muscle cells, such as fibroblasts. The cells outside 
the clusters appeared to be non-muscle cells, because 
while they exhibited active locomotion and prolifer- 
ation, but no contraction. We therefore undertook a 
histochemical analysis to determine which types of 
cells were involved in the formation of clusters. 

a-sarcomeric actin (actin expressed in muscle 
cells) and cellular nuclei were stained using a- 
sarcomeric actin antibody and DAPI (a dye that 
stains DNA), respectively. Fig. 5a shows a fluores- 
cent image of cells at the sixth day from the onset 
of incubation. Since a-sarcomeric actin is specifi- 
cally expressed in cardiac muscle cells, cell nuclei 
colocalized with a-sarcomeric actin are considered 
to be cardiac muscle cells, while those cell nuclei 
that are not associated with a-sarcomeric actin are 
not muscle cells. Therefore, using this method and 
the associated images, we are able to discriminate 
between cardiac muscle cells and non-muscle cells. 
The fraction of cardiac muscle cells, in which the 
nucleus colocalizes with a-sarcomeric actin, was 
0.76 ± 0.1. 

It was also found that the cell clusters mainly con- 
sisted of cardiac muscle cells and that non-muscle 
cells were located outside the clusters. Actually, the 
distribution of nucleus density (Fig. 5b) was simi- 
lar to that of a-sarcomeric actin. In regions where 



the local density of nuclei was large, which corre- 
sponded to clusters, the expression of a-sarcomeric 
actin was evident. From this observation, it is obvi- 
ous that the clusters consist of cardiac muscle cells. 
While there were nuclei outside clusters, these nuclei 
were not colocalized with a-sarcomeric actin, which 
indicated that these cells were non-muscle cells. The 
density of these cells was low and a region without 
cells was observed outside the clusters. 

3.3. Radial distribution of cells and the size of 
clusters 

In order to quantify the distribution of cells, the 
cells were fixed on the seventh day, and the nuclei of 
cells were stained with DAPI. The images of the cell 
nuclei were obtained using fluorescent microscopy. 
Figure 6a shows the radial distribution function of 
cell nuclei calculated from these images for several 
initial cell densities. For all densities studied here, a 
peak is observed around 20 /^m in the radial distribu- 
tion function. Considering the size of a nuclei and a 
cell (around 10 /im and a few tens /xm, respectively), 
this peak is interpreted as the distance between nu- 
clei of closely packed cells, i.e., the cells in the same 
cluster. From Fig. 6a, we considered that the nucleus 
within 50 /im from each other constitute a single 
cluster. Using this definition, the number of nucleus 
(cells) inside a single cluster was determined for ev- 
ery cluster. Figure 6b, c, and d displays the number 
distribution of nuclei in a cluster for several initial 
cell densities. As seen in these figures, the number 
in a single cluster varies from one to a thousand. As 
the initial density of cells increases, the fraction of 
large (e.g., > 10 nuclei) clusters increases. 

3.4. Effect of cell contraction 

To examine the effect of spontaneous cardiac 
muscle-cell contraction on cluster formation, we 
investigated whether inhibition of contraction af- 
fected cluster formation. For this purpose, we uti- 
lized 2,3-butanedione 2-monoxime (BDM), which 
is known to affect the actomyosin system in muscle 
fibrils and inhibit contraction [17,18]. 

Application of this agent to culture medium re- 
sulted in the inhibition of spontaneous contraction 
of cells. The activity of spontaneous cellular contrac- 
tion was estimated using the method described in 
Section 2.8. In this method, fiuorescent microbeads 
((j) 1 ^m) were dispersed on the surface of the cells 



6 




i«' IB? J# 



ig' I'H^ -# 



Fig. 6. Distribution of cell nuclei, a) Radial distribution 
function of cell nuclei for several initial cell density. The 
initial cell densities are indicated in the inset. Dotted curves 
are the fits to Eq. (2). b) - d) Histograms of number of 
cells in a single cluster. Initial cell densities axe b) 4.4 x 10^ 
cells/mm^, c) 8.8 x 10^ cells/mm^, d) 1.3 x 10^ cells/mm^, 
respectively. The average numbers of cells in a cluster are a) 
7.8 ± 23, b) 23 ± 124, and c) 40 ± 431, respectively. 

and the fraction of periodically moving beads was 
determined using video microscopy. Figure 7 depicts 
the fraction of periodically moving beads with vary- 
ing concentration of BDM. As illustrated in this fig- 
ure, the fraction of periodically moving beads de- 
creases as the concentration of BDM increases. In 
the all preparations, the presence of 10 mM of BDM 
almost completely inhibited the spontaneous con- 
traction of cardiac muscle cells. It has been con- 
firmed that the cells were alive at these concentra- 
tions of BDM, because the spontaneous contrac- 
tion resumed after removing the BDM and exchang- 
ing the culture medium with fresh medium without 
BDM (data not shown) . 

The cells were supplied with culture medium con- 
taining varying concentrations of BDM and incu- 
bated for several days. Figure 8 shows the micro- 
scopic images of cells on the fifth day of culture. This 
figure illustrates that the formation of cellular clus- 
ters is suppressed in the presence of a high concen- 
trations of BDM (> 10 mM). In order to quantify 
the distribution of cells, the cells were fixed on the 
seventh day, and the nuclei of cells were stained with 
DAPI. The images of the cell nuclei were obtained 
using fluorescent microscopy. Figure 9a shows the 
radial distribution function of cell nuclei calculated 
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Fig. 7. Effect of BDM on the spontaneous contraction of 
cardiac muscle cells, a) Moving beads assay. Fluorescent 
image of microbeads {(p 1 /im) attached to the surface of cells 
supplied with culture medium containing no BDM. Inset 

exemplifies the motion of the bead that is marked in the 
image. The scale bar represents 50 ^m. b) The fraction of 
moving fluorescent beads attached to cell surface is plotted 
against the concentration of BDM. The error bars represent 
itS.E.M. 

from these images. For cells cultivated in the absence 
of BDM, the radial distribution function exhibits a 
peak at approximately 20 /im, indicating the forma- 
tion of clusters. As the concentration of BDM in- 
creases, the intensity of the peak decreases and the 
location of the peak shifts towards a larger distance. 
At 15 mM of BDM, the radial distribution function 
exhibits almost no peaks, indicating the absence of 
clusters and a random distribution of cells. 

In order to confirm the abovementioned trends 
quantitatively, we determined the intensity and the 
location of the peaks. Since the raw data is rather 
noisy, the experimental data were fitted to the fol- 
lowing function 



g{r) = Aexp 



(2) 



7 



a) BDM OmM b) EDM 5mM 




Fig. 8. Microscopic images of cells on the fifth day from 
the onset of incubation in tlie presence of various BDM 
concentrations. The concentration of BDM is indicated on 
the top of each figure. The scale bar represents 100 /im. 

from which we determined the intensity and the lo- 
cation of the peaks. As observed in Fig. 9, this func- 
tion fits the experimental data well, and we can es- 
timate several parameters from the fitting function. 
In Fig. 9c, the intensity and the location of the peak 
are plotted as functions of BDM concentration. The 
intensity of the peak is represented as the ratio of 
the magnitude of the radial distribution function at 
the peak to the magnitude of the radial distribution 
function at a large distance R. In this case, we se- 
lected R = 200 /im. Based on these plots, it became 
clear that the intensity of the peak decreases and the 
location of the peak shifts towards a larger distance 
as the concentration of BDM is increased. 

One may notice that the average density of nuclei 
decreases as the concentration of BDM increases in 
Fig. 9a. This could be because the number of dead 
cells increases as the concentration of BDM in cul- 
ture medium increases. However, the disappearance 
of the peak in the radial distribution function with 
an increase in the BDM concentration is not simply 
attributable to a decrease in the average density of 
cells. In fact, the morphology of cells is clearly differ- 
ent for different concentration of BDM (see Fig. 8). 
Moreover, when the initial density of cells was in- 
creased three fold, a similar trend regarding was ob- 
served between the concentration of BDM and the 
shape of the radial distribution function (see Fig. 9b 
and Fig. 9c). The change of the shape of the ra- 
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Fig. 9. Radial distribution functions of cell nuclei on the sev- 
enth day. Different curves represent the radial distribution 
function at different concentrations of BDM. a) Initial den- 
sity of cells was 4.4 x lO'^ cells/mm^, b) Initial density of cells 
was 1.3 X 10^ cell/mm^, c) Intensity (top) and the location 
(bottom) of the peak radial distribution function. Experi- 
mental data displayed in a) and b) were fitted to Eq. (2), 
from which the intensity and the location of the peak was 
estimated. In a) and b), the functions, g{r), obtained by fit- 
ting are also represented using dotted curves. The intensity 
is represented as the ratio g{rinax)/g{R), with letting rmax 
be the location of the peaJj and R = 200 //.m. In the bottom 
plot, the error bars indicate the half width of g{r) — g{R). 
Crosses and squares curves represent data derived from the 
plots in a) (lower density) and b) (higher density), respec- 
tively. 
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Fig. 10. Microscopic images of cells on the substratum coated 
with varying concentrations of collagen solution. The images 
of cells fixed on the seventh day from the onset of incuba- 
tion are displayed. The scale bar represents 100 /im. The 
concentrations of collagen utilized to coat substrata are 1.2 
rag/l (a), 3.6 mg/£ (b), 120 mg/£ (c). 

dial distribution function induced by application of 
BDM is qualitatively different from a mere shift of 
the ordinate caused by the change in the cell density 
(see Fig. 6). 

Interestingly, a marked correspondence was ob- 
served between spontaneous contraction activity 
and the structure of the radial distribution function 
regarding the dependence of the contractions on the 
concentration of BDM. At a BDM concentration of 
less than 6.7 mM, both the contraction of cells and 
the formation of clusters were evident, while at a 
BDM concentration greater than 8.4 mM, both the 
contraction of cells and the formation of clusters 
were suppressed (see Figs. 7 and 9c). Based on this 
observation, it became apparent that the sponta- 
neous contraction of cardiac muscle cells was closely 
related to the formation of clusters. 

3.5. Effect of substratum 

Next, we examined the effect of the adhesive force 
of the muscle cells to the substratum on the for- 
mation of clusters. In the abovementioned experi- 
ments, substrata c;oatcd with a thin layer of colla- 
gen were used. Given that the cells were attached to 
the substratum by a collagen coating, the effect of 
the strength of cellular adhesion to the substratum 
was tested by altering the amount of collagen in the 
coating. 

We then prepared substrata coated with collagen 
at varying concentrations. The cells were plated on 
these substrata and incubated for several days under 
the same conditions as those described above. It was 
confirmed by visual inspection that the cells exhib- 
ited apparent contractions independent of the con- 
centration of collagen. Figure 10 shows microscopic 
images of cells plated on substrata coated with col- 
lagen at varying concentrations. Interestingly, clus- 
ters formed when the concentration of collagen was 



low. Conversely, no clusters formed when the con- 
centration of collagen was high. In order to quan- 
tify this observation, the radial distribution function 
was calculated. On the seventh day from the onset 
of incubation, the cells were fixed and stained with 
DAPI. The distribution of cell nuclei was then in- 
vestigated using fluorescent microscopy, and the ra- 
dial distribution function was calculated. Figure 11a 
displays the radial distribution functions of cultures 
grown on substrata with varying amounts of colla- 
gen. As seen in this figure, when the concentration of 
collagen was low (1.2 mg/^), the radial distribution 
function exhibits a clear peak around 20 /im, which 
is similar to the former experiments in the absence 
of BDM. Furthermore, the peak in the radial dis- 
tribution function became less intense and broad as 
the concentration of collagen increased. In substrata 
with a high concentration of collagen (120 mg/£), al- 
most no peak was observed in the radial distribution 




Fig. 11. a) Radial distribution function of cell nuclei on the 
seventh day. Different curves represent the radial distribution 
function at the different collagen concentrations used to coat 
the substrata, b) Parameter values obtained from fitting 
the curves in a) by Eq. (2). The functions, g{r), obtained 
by fitting are represented in a) using dotted curves. (Top) 
The intensity of the peaJj defined as g{riniix)/g{R), where 
R = 200 /im. (Bottom) The location of the peak plotted as 
a function of concentration of collagen. Error bars indicate 
the half width of g{r) - g{R). 
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function. These observations were also confirmed by 
fitting Eq. (2) to the radial distribution functions, 
which yielded estimates of the intensity and the loca- 
tion of the peaks (see Fig. lib). Consequently, while 
the formation of clusters was suppressed when the 
substratum was coated with high concentrations of 
collagen, cells were able to form clusters when col- 
lagen concentrations in the substratum were low. 

4. Discussion 

In the present study, spontaneous formation of cell 
chistcrs was observed in primary cultures of ventric- 
ular cells from newborn rats. While the formation of 
clusters has been observed in several earlier studies 
using similar preparations [19,20], the mechanisms 
underlying cluster formation have not been clarified. 
In the present study, we have obtained several new 
insights into the mechanism of cluster formation. 

First, as revealed by histochemical analysis, clus- 
ters were comprised primarily of cardiac muscle 
cells, with non-muscle cell types, such as fibroblasts, 
being located outside the clusters. 

Second, the cells in a cluster are almost fully 
packed. The number of cells incorporated in the 
same cluster increases as the density of cells in- 
creases. 

Third, the direction in which the cells move is cor- 
related to the direction of the long axis of cells and 
the direction of cell contraction (Fig. 4). Further- 
more, the observation that the inhibition of spon- 
taneous contraction of cardiac muscle cells resulted 
in suppression of cluster formation (Figs. 8 and 9) 
supports the idea that the spontaneous contraction 
of cardiac muscle cells is related to cluster forma- 
tion. It is thus plausible that the periodic contrac- 
tion generated tension between the cells, and that 
this tension makes the cells move. 

Fourth, the concentration of the collagen coated 
on the substratum affected cluster formation. When 
the substratum was coated with a solution contain- 
ing a high concentration of collagen, formation of 
cell clusters was suppressed (Figs. 10 and 11). The 
adhesive force of cells to the substratum is expected 
to increase when the concentration of collagen used 
to coat the substrate increases, because the number 
of contacts formed between the cytoskeleton and col- 
lagen would increase. Consequently, cluster forma- 
tion is suppressed when the adhesive force between 
cells and the substratum is large. 

Based on the above observations, the mechanism 



of cluster formation by cardiac muscle cells is pro- 
posed, as schematically shown in Fig. 12. In the first 
stage, the distribution of cells is random. At the ini- 
tial cell density (4.4 x 10^ cells/mm^), the area of 
substratum covered by cells is approximately half. 
Consequently, individual cells are likely to being in 
contact with several neighboring cells. Once two car- 
diac muscle cells establish contact, gap junctions are 
formed [11], and their contractions become synchro- 
nized (see Fig. 2). These cells are also mechanically 
coupled via intercellular contact, including adher- 
ence junctions and dcsmosomes [22]. Contraction of 
a cell is attributed to muscle fibrils, which are com- 
plexes of cytoskcletons. The contraction of the mus- 
cle fibrils exerts tension on cytoskcletons, including 
actin filamc^uts, which are anchored to the substra- 
tum via molecular apparatuses termed focal con- 
tact (see Fig. 12a). Owing to this tension, focal con- 
tacts at the periphery of cells are pulled into the 
center of each cell (this process may involve rapid 
destructions and reconstructions of focal contacts). 
Through this process, cytoskcletons are rearranged 
and become compact. Intercellular contacts are also 
expected to play an important role in maintain- 
ing and stabilizing the mechanical contacts between 
cells, against the contractile tension of cytoskcle- 
tons [21]. If these intercellular connections are suf- 
ficiently tight, the center of each cell become closer 
to each other because each cell becomes compact 
(see Fig. 12b). Since a number of cells are mechani- 
cally connected to each other and each cell becomes 
compact, the population of cells becomes denser and 
forms a cluster. 

If the force generated by contraction is too weak 
to move focal contacts, contraction of cells will not 
result in the cells moving closer. Similarly, if the con- 
centration of collagen is large, the number and the 
strength of focal contacts between the cytoskeleton 
and collagen would increase, making rearrangement 
of the cytoskeleton less likely and leading to the sup- 
pression of cluster formation. Furthermore, as the 
density of cells increases, the number of cells that are 
initially in contact increases, resulting in the larger 
size of clusters. This scenario is consistent with the 
experimental results on the effects of BDM and also 
of the collagen concentration. 

With regard to the cluster formation in cardiac 
cell culture, it has been reported previously that 
the presence of fibroblasts is important for induc- 
ing the aggregation of cardiac muscle cells in some 
situations [21]. In the present study, given the diffi- 
culties associated with complete removal of prolifer- 
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Fig. 12. Hypothetical mechanism of cluster formation, a) 
Schematic cross section of a cardiac muscle cell. Cells contain 
muscle fibrils, which become connected to collagen via focal 
contacts at their ends. The substratum is coated with colla- 
gen. Contractile motion is indicated by arrows and adjacent 
cells are connected through intercellular contacts, including 
adherence junctions and dcsmosomes. b) Dorsal view of two 
cardiac muscle cells in contact with each other via an inter- 
cellular contact. Individual cells contain muscle fibrils and 
exhibit periodic contractions, resulting in the focal contacts 
at the ends of the muscle fibrils experiencing periodic force. 
These actions lead to the rearrangement of focal contacts 
and cytoskeletons in individual cells, which brings the two 
cells closer together. 

ating non-muscle cells from cultures, it is not clear 
whether the presence of non-muscle cells, whose 
fraction was not large in the present study (around 
20 % even in the late stage of culture), is essential 
for the observed phenomena. However, as discussed 
above, several observations in the present study 
suggest that formation of clusters is also mediated 
by the interactions among cardiac muscle cells. 

Finally, let us discuss the relation of the present 
results to the previous studies on the wave dynamics 
on the cardiac tissue culture. The formation of cel- 
lular clusters causes the cell preparations to become 
inlioniogcnoous and this may alter the dynamics of 
excitation waves [3]. In particular, Steinberg et al. 
reported that the propagation of excitation waves 
in a cardiac cell culture is affected by the concen- 
tration of collagen used to coat the substratum [23]. 
Since we have observed that the formation of clus- 
ters is also affected by the concentration of collagen 
on the substratum, it is possible that the alternation 
of wave dynamics observed in Ref. [23] is directly re- 
lated to cluster formation. Furthermore, while it has 
been demonstrated previously that the dynamics of 
excitation waves vary in a complicated manner over 
the time scale of several days [7,8,9,10], the results 



presented here demonstrate that that formation of 

clusters is one of the essential factor to cause these 
differences in observed wave dynamics. 

5. Concluding Remarks 

In the present paper, the mechanism of sponta- 
neous cluster formation in a primary culture of ven- 
tricular cells from neonatal rats was investigated. 
On the basis of several observations, it is suggested 
that periodic contraction of cardiac muscle cells re- 
sults in the relocation of focal contacts, and that this 
causes cells to aggregate. In other words, the evolu- 
tion of a cellular network is directly affected by the 
contraction activity of individual cells. 

The mechanism described here is still hypotheti- 
cal, and further elucidation is necessary. By visual- 
izing cytoskeletons and focal contacts in living cells 
using fluorescent proteins, more detailed informa- 
tion on the molecular processes taking place during 
cluster formation would be obtained. Construction 
of theoretical models would also facilitate our un- 
derstanding of the molecular mechanisms involved. 

In addition to the molecular mechanism under- 
lying cluster formation, it would be interesting to 
investigate how the formation of cluster structure 
affects the contraction dynamics of cells. Moreover, 
the effect of spontaneous cell contraction on cellu- 
lar properties, such as shape and the orientation of 
cells should also be investigated. 

In previous studies, the contraction dynamics of 
cardiac cell cultures have been theoretically stud- 
ied within the framework of conventional reaction- 
diffusion systems or coupled oscillators [3,5,7,9,10], 
with age-dependent changes in the cellular network 
considered as merely changes in external parame- 
ters. The present study, however, suggests that the 
changes in these parameters may in fact be induced 
by the dynamics of the contraction itself. It might 
therefore be possible that a new class of phenom- 
ena arises through interactions between the dynam- 
ics of contraction and the structure of the cellular 
network. This research may therefore lead to the de- 
velopment of new strategies for controlling the ar- 
chitecture of cellular networks by modulating the 
contraction activity of cells. 
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